Influence of the argon admixture on the reactive oxide species formation inside an atmospheric pressure oxygen plasma jet

In this work, a new atmospheric pressure plasma generated in a wire-to-multiwire dielectric barrier discharge on pure oxygen is introduced. This special geometry of 13 wires (one central wire and 12 ones on the external tube) is feeding by a radio frequency (RF) power (13.56 MHz, 1 kW) and produces a stable discharge. The capacity of this device to produce oxygen reactive species and the influence of Ar gas mixture (1–3%) on this production are investigated. The main characteristics of this DBD plasma are measured using optical emission spectroscopy techniques. The rotational, vibrational, and excitation temperatures along with the electron density are determined from OH (A2Σ → X2Π) band and the Stark broadening of the hydrogen atomic line at 486.1 nm, respectively. The temporal evolution and spatial distribution of charged and reactive species in this plasma are also numerically studied by a Global scheme and a two-dimension fluid model based on drift–diffusion approximation. A kinetic dominated by electron collisions is obtained for this plasma. The generation and movement of electrons, positive and negative ions in the wire-to-multiwire configuration are analyzed and discussed according to changes the electric field and plasma frequency. It is shown that the density of both charged and reactive species increases by adding a small amount of argon to the oxygen plasma while the electron temperature reduces in this configuration. A high level of agreement is observed between the experimental and simulation results for the electron density and temperature in this DBD plasma.

It is well-known the high number of applications of the oxygen plasma, from the treatment and cleaning of surfaces to the biomedical field.The various components of plasma, such as electrons, ions, neutrals, molecular and atomic species, interact with the surface, altering its chemical and physical properties.Oxygen discharges have been applied for changing the properties of the surface of metals, plastics, glass or polymers 1,2 .Polymers are becoming increasingly popular due to their affordability and superior performance.Because of their low surface energy, which results in poor adhesion, their applications are limited.Functional groups formed by the oxygen plasma, are introduced into the surface, producing its modifications without disturbing the properties of the material 3 .
Surface modification can significantly enhance the adhesion of the polymer surface.Moreover, prior to growing an oxide mask layer or depositing an antireflection coating, it is essential to ensure that the surface is thoroughly cleaned in the case of microelectromechanical systems.The hydrophilic/hydrophobic properties of the surface or its poor adhesivity can be improved by the plasma treatment 4,5 .An atmospheric pressure oxygen plasma could be a convenient, in-line treatment process for cleaning, modifying, and activating the surface.So, it has been applied for photoresist stripping, removing polymer films, and oxidation or deposition of oxides film 6 .
Plasma treatment is a highly effective method for cleaning contaminants and modifying surfaces, making it ideal for a variety of applications.The oxygen plasma is also used for sterilization and cleaning of surfaces.The OPEN 1 Faculty of Physics, Shahid Bahonar University of Kerman, Kerman, Iran. 2 Department of Physics, Vali-e-Asr University of Rafsanjan, Rafsanjan, Iran. 3 Photonics Institute, Kerman Graduate University of Technology, Kerman, Iran. 4 Faculty of Civil Engineering and Environmental Sciences, Bialystok University of Technology, Bialystok, Poland. 5Department of Physics, School of Engineering Sciences of Belmez, University of Cordoba, Cordoba, Spain.* email: alibarkhordari20@yahoo.com

Experimental setup
As depicted in Fig. 1a for the designed RF plasma discharge, a central tungsten wire (inner) with 4 mm diameter has been placed at the center of a quartz tube with inner and outer diameters of 10 mm and 12 mm, respectively.Spiral-shaped roughness has been made in the inner wall of the quartz tube to move the gas as a vortex flow in the discharge medium.Due to this movement of gas flow, the temperature of the discharge medium is decreased by this roughness dielectric barrier.Moreover, 12 aluminum wires of 1 mm diameter have been cylindrically fixed around the tube as ground electrodes, with a separation of 30°.This device operates at atmospheric pressure with pure oxygen and small argon amounts mixed by oxygen.Thus, the plasma discharge is formed in a 1 mm gap between tungsten inner rode and aluminum outer wires that are covered by quartz dielectric.The length of generated RF plasma is about 20 mm along the inner electrode.This arrangement for plasma jet strongly prevents arcing between the inner and outer electrodes which is particularly unwanted for biomedical applications.Such a configuration for plasma jet generates a stable plasma discharge.The produced plasma transfers the reactive species to the processing region, which avoids disturbance to plasma stability.
In this work, a power supply with RF fixed-frequency (13.56 MHz) and constant output power (1 kW) has been applied.At first, both the output RF power and impedance are varied to reach the most stable plasma discharge for pure oxygen and all the its mixtures with argon.The impedance matching network works based on a transmission line and, a lumped L-type network has been used.It has to be noted that, the L-type matching network is lossless or, at least, the loss can be made extremely small with the proper component choices.However, the RF input power depends on the application and dimensions of the RF plasma jet.
At a DBD (AC or pulsed) plasma with high voltages, a breakdown occurs if the maximum voltage is reached.At atmospheric pressure, the streamer breakdown is contracted.Then, thin discharge channels appear, called filamentary.At RF discharges, penning ionization, stepwise ionization, and charge transfer collisions have a crucial effect on ionization dynamics that avoid the formation of this filamentary.So, a diffusive DBD or glow DBD is formed 19,20 .
In addition, two spectrometers are used for the spectroscopic characterization of this discharge.A HR4000 Ocean Optic spectrometer with optical resolution of width ~ 0.025 nm was applied for spectroscopic measurements.It stores a full spectrum at every millisecond with a wavelength range sensitivity of 200-1100 nm.Also, a spectrometer with higher resolution, Jobin Yvon, TRIAX550 with a resolution of 0.05 nm at 500 nm is used for the determination of plasma parameters, electron density and temperatures.It must be noted that, in order to calibrate these spectrometers, a HG-1 Mercury-Argon lamp was used as a light source that is capable of producing the spectral emission lines.Moreover, an optical fiber probe was placed at the middle of the plasma region (10 cm) which was perpendicular to the plasma discharge with a spatial distance of 5 mm.
It should be mentioned that the operational parameters of the external circuit are fixed.Using the gas mixer, the pure oxygen with 15 SLM (Standard Liters per Minute) gas flow rate, and O 2 /Ar gas mixture by 15 and 1-3 SLM gas flow rates enter the plasma jet during the ignition.When stable discharge is produced, a lower working flow of 1 SLM is used.It is worthy to mention that this flow rates are controlled using a velocity gas flow controller (Fig. 1b).

Modeling
In this work, the effects of various mixtures of oxygen and argon on the spatial and temporal distributions of charged and reactive species, neutral particles, and the electron temperature in the wire-to-multiwire dielectric barrier discharge are theoretically studied when an RF power is used.The global scheme and 2D model are respectively based on the kinetic and fluid approaches which have been used to simulate O 2 /Ar gas mixture in the introduced electrode configuration.

Global model
Let us consider a cylindrical chamber with a radius of R and a length of L into which the neutral species is entered by a steady flow via the gas inlet.It is assumed that the contents of the chamber are uniformly distributed at space, and the energy is uniformly deposited into the plasma.Moreover, electrons possess a Maxwellian-like energy distribution.Investigations of the electron energy distribution function (EEDF) reveal that it is largely consistent with the Maxwellian distribution in molecular plasmas 21 .
A continuity equation is used to describe the creation, gas phase and surface reactions, and losses of each species.So, a system of first order differential equations should be solved to evaluate the formation of reactive oxygen species originating from the oxygen molecules and the influence of adding argon atoms on the plasma processes in wire-to-multiwire dielectric barrier discharge considered in this work.The quasi-neutrality condition requires that the charged particle species be in balance: where n refer to the number density of electrons and different positive and negative ions of oxygen and argon species considered in this study (see Table 1).In general, the time evolution of the number density of different species, n s , is expressed as [22][23][24] : (1) where N s is the number of reactions influencing on s specie, R i,s and R s,j refer to the reactions rate of the reactions that populate and depopulate this specie.Moreover, the reaction rate R i,j between i and j species, depending on the density of each species, n i and n j , with k i,j , the reaction constant, as follows 25,26 : In this work, the reaction rate constants for electron collisions are calculated using the collision cross sections, σ i,j (ε), and EEDF, f(ε), by the following integral 27 : with ε and m e being the energy and electron mass.Otherwise, the reaction rate constants for heavy particles collisions are calculated by: with T g being the gas temperature.It is necessary to note the neutral species are diffused in the considered cylindrical chamber, and their diffusion coefficient is expressed as 28 : where T g is the gas temperature, m n is the mass of the neutral species, v n = (8eT g /πm n ) 1/2 is the mean speed of the neutral species, and i is the mean free path that is given by 29 ; with n g,j and σ i,j being the neutral species of jth ion and the scattering cross section of ion-neutral pairs for jth neutral species.In addition, the effective length of the diffusion of each neutral species in the cylindrical chamber with a length of R and a radius of R can be defined as follows 28 : These quantities allow us to find the effective loss-rate coefficient required for examining the diffusional losses of the neutral species such as oxygen atoms, excited argon atoms, metastable oxygen atoms and molecules that reached the jet wall 29 ; where V denote the volume of the cylindrical chamber, A is the surface area of the wall, and γ n being a coefficient for the neutral species stuck on the wall surface.On the other hand, the rate coefficient of the ion flux could be described as 29 : where u B = √ eT e /m i refers to the Bohm velocity, with m i being the ion mass.Furthermore, the quantity of h L and h R denote the ratios of the positive ion density at the edge to center of the cylindrical chamber of plasma respectively given by 29 : (3) www.nature.com/scientificreports/with α = n − /n e is the electronegativity with n − being the negative ion density, and γ = T i /T e .For the regime of γ (R, L) ≤ i ≤ (R, L) , the above equations tend to Godyak's equations for an electropositive discharge because of α → 0 30 .In order to stabilize the modeling, the electron energy density equation should be solved together with the number density equation in simulation procedure which is described by 31 : where P is the input power density, Q elas and Q inelas are the energy loss and net energy loss due to the elastic and inelastic processes, respectively 32 .The power balance equation calculated by the ratio of the absorbed power P abs and the power losses resulting from elastic and inelastic collisions as well as the charged particle flow to the walls, is defined as 33 : c being the energy loss per electron-ion pair produced by the neutral which is given as 33 : where T e is the electron temperature, m e is the electron mass, m i is the ion mass.The coefficients of k el , k iz , and k ex,i are the rate constant of elastic scattering, ionization, and ith excited state, respectively.Moreover, E ex,i and E iz refer to the energy of the ionization and the ith excitation processes, respectively 33 .
The system of these equations was solved Global Module of PLASIMO.It should be mentioned that this module works with BOLSIG + to use of the needed data of collisional electron reactions 32 .This model gives the time evolution of number density of k species that can be as a criterion of the lifetime of generated species in the plasma discharge.
In this model, 56 species presented in Table 1 have been considered which are resulted in 4000 reactions including electron impact, electron attachment, electron-ion recombination reactions, neutral-neutral reactions, ion-heavy particle reactions, vibrational energy transfer reactions.

Model domain
Figure 2 shows a schematic of the domain which has been considered for simulation.As seen, a central tungsten wire with a diameter of 4 mm placed at the center of a 1 mm thick quartz tube was used as the power electrode.Moreover, 12 aluminum wires (2 mm diameter) were cylindrically fixed around the tube as ground electrodes.It is supposed that this structure operated at atmospheric pressure using both pure oxygen and a mixture of oxygen with small amounts of argon.Thus, the plasma formed in a gap of 3 mm between the tungsten inner www.nature.com/scientificreports/rode and the aluminium outer wires covered by the quartz dielectric.Practically, this arrangement for discharge efficiently prevents arcing between the inner and outer electrodes, which is particularly unwanted for biomedical applications.Such a jet-type configuration generates a stable plasma.The produced plasma transfers the reactive species to the processing region, which avoids disturbance to plasma stability.Due to cylindrical symmetry in the presented structure, the cylindrical coordinate is used to spatially simulate the 2D fluid model.To this aim, a symmetry along the z-direction is considered by assuming being in the middle of discharge (a large distance from two end of wires).Then the spatial description of the problem is only defined in radial and azimuthal directions.Because of using the fluid model to spatially simulate, only three charged species (electrons, positive and negative ions) are taken into account.In following, the fluid and the global models are briefly described.

Model equations
The spatial and temporal macroscopic description of the gas discharge inside the wire-to-multiwire configuration is via solving the fluid continuity equations for different species coupled with Poisson's equation.These equations based on drift-diffusion approximation are solved using the finite difference (FD) method.The flux for each species in the drift-diffusion approximation is based on the momentum conservation of each species which is defined as 29 : where μ i is the electrical mobility, D i is spatial diffusivity, E is the electric field, and n i is the number density of species i.The plus or minus sign in this relation accounts for the sign of the charged particles 34 .The continuity equation for all the formed species in the plasma discharge is expressed as 25 : where n i is the number density, Γ i expresses the flux for species i and R i,m is the reaction rate between species i and species m.In this model, the kinetic scheme is the same that the used for global model (see table in supplementary material).
Moreover, the rate of change of the electron energy density is described by 29 : where n ε is the electron energy density, eE • Ŵ e is the ohmic or joule heating for electrons, R ε is the energy loss or gain due to inelastic collisions which is obtained by summing the collisional energy loss or gain over all reactions 25 : where x j is the mole fraction of the target species for reaction j, k j is the rate coefficient for reaction j, N n is the total neutral number density and Δε j is the energy loss from reaction j.
For non-electron species, the following equation is solved for the mass fraction of each species 35 : where ρ denotes the density of the mixture, ξ k is the mass fraction of the kth species, u is the mass averaged fluid velocity vector, V k and R k are the multicomponent diffusion velocity and the rate expression for species k, respectively 35 .In a diffusion model, the multicomponent diffusion velocity, V k , can be calculated by the Fick's law 35 : where D k,f and D T k are the diffusion coefficient and thermal diffusion coefficient for each species, M is the mean molar mass of the mixture, T is the gas temperature, and z k and µ k,m represent the charge and mobility of species k.This mobility is given by the Einstein's relation 35 : The energy balance for heavy particles is calculated by Fourier`s equation that allows to obtain gas temperature T: where k and C p are the thermal conductivity and specific heat capacity of oxygen gas, respectively, and R is the heat sources, originating from the Joule effect and collisions among electron and heavy particles 35 .In the model, ( 16) the axial energy transport is considered to be neglected because low working gas flow is used.It should be highlighted that in a discharge of molecular gas, the fast-heating effect by quenching of excited molecules with neutral or atomic species could be vital in plasma heating.It becomes more crucial in air discharge by quenching N 2 excited states and oxygen molecules [36][37][38] .This effect is neglected in the studied plasma discharge because the amount of nitrogen is too weak.Although the quenching of O( 1 D) oxygen occurs in pure oxygen plasma, this effect is relevant for high energy density, as in the case of nanosecond discharge 37 .The RF discharge is far from these conditions.So, we can neglect the fast heating in our model.
To initiate dielectric barrier discharge in the wire-to-multiwire structure, electric potential should be applied between the electrodes, thus the Poisson's equation must also be considered in the model 25 : where φ is the electric potential, ε(= ε 0 ε r ) is the permittivity ( ε 0 for vacuum and ε r for the dielectric), n k and Z k e denote the density of charged species and their charge, respectively 39 .

Boundary conditions
To obtain a unique solution for the coupled equations system with the geometry presented in Fig. 1, necessary boundary conditions must be imposed.The applied boundary conditions for the wire-to-multiwire dielectric discharge are similar to those that can be found in the existing literature 40 .The following boundary condition is used to account for the particles flux in wall: where n is the normal vector pointing toward the tube wall and, v th(i) is the thermal velocity of particles 41 : and the number a i is defined by: For electrons, as a special case, the particles flux due to secondary electron emission (SEE) is added to the system and is defined as follows 40 : where γ p are the SEE coefficients, which defines the average number of electrons emitted per impact of ions p on the tube wall.Similarly, boundary condition for electron energy is 41 : Here, the second term is the SEE energy flux, being ε p the mean energy of the secondary electrons.The plasma discharge is driven by an RF potential applied to the wire centered at the axis of quartz tube as and other electrodes located on the tube are grounded.It should be mentioned that ϕ RF is the amplitude of RF voltage used to initiate the plasma discharge.The electric displacement vector, D = εE , on the dielectric tube is changed as follows: and the surface charge density, ρ s , on the dielectric surface is given by: where j e and j i are the total electron and ion current densities on the dielectric layer 40 .Unless otherwise mentioned, the simulations have been performed with the chosen parameters presented in Table 2.It has been tried to select the parameters close to the real experimental values of atmospheric plasma discharge.The reactions used in the 2D fluid model have been presented in Appendix (see Table S1).It must be noted that the drift-diffusion module of PLASIMO has been applied to do the 2D simulation.

Spectroscopic results
By analyzing the spectra recorded at various percentages of Ar in the O 2 /Ar gas mixture, different species generated in the wire-to-multiwire DBD plasma at atmospheric pressure are determined.Figure 3 shows both typical emission spectrum of pure O 2 and O 2 + 3%Ar plasmas and 1 kW RF power.It must be noted that the tube of this plasma reactor was opened at its end, which allowed some air went into it and reached the plasma.Thus, some excited nitrogen-containing species were detected in the spectra including N 2 and NO molecules, N 2 + ions, and nitrogen atoms.This nitrogen entrance could produce heating of the discharge, especially in the outside flame, and influence in the application of this type of plasma.We are considering neglected this effect because we study the oxidative specie formations inside the reactor.Water molecules are usually present in ambient air and in the gas as impurity, which explains the presence of OH radicals and hydrogen atoms in the spectra.On the other hand, some excited tungsten and sodium atoms were also detected, revealing that some tube and inner electrode etching took place.The gas temperature in this reactor is near to the ambient temperature (300 K).Then this etching could only be explained by the collisions of accelerated electron and ions with the quartz tube.Details of the different excited species produced in the plasmas are gathered in Table 3.Moreover, variations of the intensity of O 2 (Schumman-Runge), O 2 (2nd Neg), O I, and Ar I species at the wavelength of 266.33 nm, 254.55 nm, 777.19 nm, and 763.51 nm are represented in Fig. 4. The (H α ) and (H β ) visible spectral lines in the Balmer series of hydrogen atoms have been detected at the wavelengths of 656.1 nm and 486.01 nm, respectively 42 .These lines can be used to determine the electron density in the plasma discharge by calculating Stark broadening which is following discussed in details.

Rotational, vibrational and excitation temperatures
The electron and ions are responsible of energy transferring from the electric field.Electrons get the energy from this field and transferred to neutrals and ions by elastic and inelastic collisions.Inelastic collisions with the www.nature.com/scientificreports/molecules can produce the excitations of the ro-vibrational states, and also the excitation of electronic levels of these molecules if the electrons have enough energy.The rotational, vibrational, and excitation (or electronic) temperatures, T rot , T vib , and T exc , give the distribution of rotational, vibrational and electronic states of the molecules, considering a Boltzmann distribution.These temperatures can be obtained using the OES technique which is characterized from the energy levels of excited states.To examine the rotation, vibration, and excitation temperatures in the wire-to-multiwire DBD plasma, the OH ( A 2 → X 2 ) band is chosen in the emitted spectra of the plasma discharge at atmospheric pressure.Hence, the simulated spectra are obtained and, to estimate the rotation, vibration, and excitation temperatures, they are fitted with the experimental emission spectra using the SPECAIR software 43,44 .In this simulation work, all the affected parameters are considered as the instrumental resolution and line shape with the collisional broadening.
Using the ro-vibrational spectrum of •OH radicals at the wavelength of 309.1 nm ( A 2 → X 2 transition), the rotational, vibrational, and excitation temperatures can be obtained 45 .Variations of all temperatures in terms of Ar contributions in the O 2 /Ar gas mixture are illustrated in Fig. 5a.As can be seen, all temperatures are intensively reduced by adding Ar to the O 2 plasma and these reductions are continued by increasing Ar amounts in the O 2 /Ar gas mixture.The values of T rot and T vib are lower than the T exc .In equilibrium conditions the T rot is Table 3. Detected peaks for the RF plasma jet with O 2 /Ar gas mixture 42,43 .www.nature.com/scientificreports/assumed close to the gas temperature T g and the T exc close to the electron temperature, T e ,.In non-equilibrium conditions, the following relation is followed:

Electron density
The electron density of this plasma was determined from the OES technique based on the Stark broadening of hydrogen Balmer series H β line measurement.For this determination, they considered the different interaction mechanisms causing the broadening of this line, whose profile can be considered Voigt-type, resulting from the convolution of a Gaussian and a Lorentzian profile 46 , where M is the atomic weight of hydrogen in atomic mass unit, T g the gas temperature in K, and λ = 486.1 nm.The instrumental broadening, Δλ I , due to the spectroscopy detection system was considered equal to 0.05 nm.
Both determine the broadening of the Gaussian profile of the line whose broadening is given by 47 : On the other hand, the collisional mechanisms leading to the Lorentzian profile in this case are: (i) the Van der Waals broadening, Δλ van der Waals , given by the following expression calculated as in ref. 41 , considering for this case H 2 O molecules as perturbers of H atoms emitters 42,46 : where µ is the reduced mass of the colliding particles (µ = 0.94737 for H/H 2 O collision), λ is the wavelength of the H β line (486.1 nm), α is the molecular polarizability of the H 2 O disturbing particles (α = 1.43 × 10 -24 cm 3 ), 〈R 2 〉 is the difference of the square radius of the upper and lower levels of H β transition, T g is the gas temperature in K, and P is the pressure (1 atm for atmospheric pressure).
(ii) the Stark broadening, Δλ Stark , which is related to the electron density 48 : The resonance broadening, due to the dipole-dipole interactions of the emitters with the ground-state atoms of the same element, can be neglected in this particular case where just small amounts of hydrogen are present in the plasma.Thus, the Stark and van der Waals broadenings determine a Lorentzian profile whose broadening is expressed as 47 : Finally, every emission line profile is considered as Voigt shaped (resulting from the convolution of the Gaussian and Lorentzian functions) with a broadening given by 48 : www.nature.com/scientificreports/Therefore, a numerical fitting of the experimental spectrum of the line to a Voigt function allows us to discriminate Lorentzian and Gaussian contributions.Thus, once measured Δλ L and determined Δλ VdW contribution from (36) (using the value of T g previously measured), the Stark broadening of H β , and consequently the electron density, can be obtained (Eq.37).In this work, the Microcal Origin software was used to perform the fitting of the line to a Voigt profile.As shown in Fig. 5b, the electron number density increases at the higher argon contributions.In O 2 /Ar gas mixture, part of the energy applied to the discharge is used for the excitation and decomposition process of O 2 molecules.Thus, ionization and electron production are reduced when the concentration of these molecules is increased.

Formation of species in the plasma
The temporal evolution and formation of species in the plasma discharge was studied by the Global Model. Figure 6a shows the results of the temporal variations of the electron density and temperature for pure oxygen plasma and for different Ar admixture.As can be observed, the electrons start to be produced in earlier times when the Ar is added.It is due to the argon needs less electrical energy for its ionization.When the electrons are produced, they accelerate by the electric field and produce new ionizations.By the same motive, at higher contents of argon in the O 2 /Ar gas mixture, the number density of electrons produced in the plasma discharge is increased.In addition, the electron number density evolution presents a maximum value in 0.1 ms and then drops to reach a stable value in all considered mixtures.The increase of electron number density in the plasma discharge results in more dissociation and creation of different oxygen species due to collisional ionization processes.This behavior and density values agree with the experimental measures of electron density obtained in Fig. 5b.
On the contrary, the electrons temperature is significantly reduced at higher fractions of argon in the O 2 / Ar gas mixture as shown in Fig. 6b.It is due to the energy supplied by the electric field is distributed between the higher number of electrons.Like the electron density, the temporal evolution of electron temperature has a maximum in earlier times.But the position of this maximum is in shorter times than the electron density case.When the electron number starts to increase, the average shared energy is lower and consequently the electron temperature decreases.Also, the temperature starts to rise in earlier times as the Ar concentration is increased.The behavior of electron temperature agrees also with experimental measurements.The temperatures decrease when Ar is added to the discharge.However, values of this temperature are higher than those obtained ones by spectroscopic measurements (see Fig. 5a).This difference is due to the separation of the excitation states from the electronic states of molecules by electrons at non-equilibrium conditions, shown in Eq. (33).
Figure 7 shows the temporal changes of the positive ions number density, i.e., O + and O + 2 , at various contributions of argon added to oxygen plasma.At higher percentages of argon in the O 2 /Ar gas mixture, the number density of O + 2 is decreased while the most amount of the O + number density is resulted in 2% of Ar in the gas mixture.Besides, the number density of O + species has been grown up at all gas mixtures in respect to the pure oxygen.It must be noted that, in addition to increasing of the O + and O + 2 production, the permanence time of these species is significantly increased by adding argon gas to the oxygen plasma.As can be seen, the durability is also raised at more contents of argon in the gas mixtures.
This temporal behavior is consequence that the main sources of the O + and O + 2 positive ions production are the collisions of oxygen molecules with electrons (Reactions 11 and 13 in Table S1): So, this temporal distribution of formed ions has a direct relation with the number of electrons in the discharge.Nevertheless, the position of the maximum is a bit lower than the electron number density.This is due to the recombination and dissociative recombination processes by electron collisions start to be also important when the electron density is enough high (Reactions 20 and 48 in Table S1) and the positive ions density is consequently reduced: The temporal changes of the number density of the negative ions of oxygen atom and molecule, i.e., O − and O − 2 at different argon contents added to oxygen plasma are shown in Fig. 8.When argon is added to the pure oxygen, a broader peak is again obtained and the O − and O − 2 species are generated earlier.However, the maximum number density of the O − and O − 2 species is lightly lower at the O 2 /Ar gas mixtures than the pure oxygen and the maximum position is earlier than the positive ions (see Fig. 7).
Again, the negative ions O − are mainly produced by electrons collision with oxygen molecules (Reactions 8 and 9 in Table S1): www.nature.com/scientificreports/Then these ions start to be generated as the electron density is increased.But their reactions with: Electrons: Positive ions: And neutral species: can produce their losses.Consequently, the density of negative ions decreases when the amount of all species is enough high.Also, Fig. 8 shows that the shape of the temporal evolution of positive ions is more influenced by the Ar concentration than previous cases.It is due to the negative ions can be also destroyed by collisions with ions Ar + : The O − 2 species are generated by the collision between the negative ion of O atoms and excited species of O 2 molecules in addition to the electron collisions with the neutral species of O 2 molecules as: The values of O number density at different contents of argon added to oxygen plasma are represented in Fig. 9a.The O is produced mainly by the dissociation of the oxygen molecules by electron collision: and recombination of the O − negative ions, due also to the electron collisions: www.nature.com/scientificreports/So, it presents a maximum as the electron starts to increase and decrease similarly to the electron density.It can be seen that as soon as the oxygen atoms are created, their positive ions begin to produce in the plasma discharge by ionization and as they increase, the oxygen atoms are reduced: Again, this maximum moves to the lower times when the Ar concentration in the O 2 /Ar mixture is increased.Figure 9b shows the temporal variations of number density of neutral reactive species, O and O 3 .It should be noted the high ozone production in this plasma discharge.This specie presents the highest density of all oxidative species formed with a maximum value near 10 10 cm −3 .This fact is very interesting for the possible application of this plasma discharge.As can be seen, the number density of O 3 species has lower peak when the argon is added to the oxygen plasma, in respect to the pure oxygen.But, the durability of the ozone species is significantly increased at more contents of argon in the O 2 /Ar gas mixture.Also, the ozone species in the oxygen plasma mixed by argon gas is produced very earlier than the pure one.The ozone species is produced by the atomic collisions of oxygen atoms and ions with molecules: and destroyed mainly by the electron dissociation:

Spatial distribution of species in the plasma
The study of the spatial distribution of plasma parameter and species allow us to know where the species are formatted and to understand why this reactor design is suitable for the production of reactive oxides species.A 2-D distribution of species and plasma parameters are calculated for times 1 µs, where the plasma starts its formation and 1 ms, and we can consider that is it stabilized.
Figure 10 shows the variations of electrons density in the wire-to-multiwire DBD plasma with pure oxygen after 1 μs and 1 ms.The number density of electrons has a sharp peak close to the power electrode at the time of 1 μs and then, the number density drops in this region and the electrons travel along the radius to the center of tube as observed in the time of 1 ms.This is due to the fact that the electron collision ionization and dissociative ionization impacts are the main processes that create the electrons in this DBD plasma.Moreover, the crosssection of the collision and dissociative of electrons is more related to the electron temperature.The profile of electron temperature has a maximum value close to the ground electrode covered by dielectric in respect to other points near to dielectric (see Fig. 12), so the number density of electrons is higher in this region than other points on the dielectric.It should be mentioned that the electrons are easily diffused in the whole tube volume due to their lighter weight in respect to other species with the diffusion coefficient D e = k B T e /m e v c , where k B is the Boltzmann constant, T e is the electron temperature, m e is the electron mass, and v e is the momentum transfer collision frequency (see Fig. 13).
Figures 11 represents the spatial variations of the electron number density in the plasma discharge medium at different fractions of argon added to the oxygen plasma.It is shown that the electrons number density is increased at more contributions of argon in the O 2 /Ar gas mixture and it reaches to a maximum in a ring shape near the power electrode while it drops at the closet to the power electrode.The argon atoms in the meta-stable state, Ar * , is one of the important species to additionally produce electrons in this case of plasma that is quenched by O 2 , and hence the electrons number density is increased by more adding Ar to the O 2 plasma.On the other hand, . Spatial variations of density profiles for electrons at 1 μs and 1 ms with pure oxygen.
the attachment process and dissociative attachment reactions with neutral radicals decrease the electrons in the plasma with pure oxygen gas and smaller contents of Ar in the O 2 /Ar gas mixture.Moreover, the profile of electrons temperature is shown in Fig. 12 at 1 kW RF power with pure oxygen and various percentages of argon mixed by oxygen.As can be seen, the electrons temperature is notably decreased by adding argon to the oxygen plasma.The electrons temperature has a peak close to the power electrode and it gradually drops along the radius of tube except close to the ground electrodes which is higher than other points on the dielectric.The electric field induced by RF power is stronger close to the power electrode (see Fig. 18), and thus the electrons get more energy there, and this results in the more electron temperature in this region.However, the electron temperature is decreased in the regions which the electron collision frequency is higher and therefore the diffusion is lower, because of the lack of electron heating.Besides, the loss of collision energy per electron-ion pair reduces with adding and increasing the argon gas to the oxygen plasma discharge.Furthermore, the electrons will receive the additional energy loss at argon contributions in the O 2 /Ar gas mixture, and so the electron temperature grows up.
Figure 13 illustrates the collision frequency of electrons in the wire-to-multiwire dielectric barrier discharge at the RF power of 1 kW with pure oxygen and different contents of argon gas added to the oxygen plasma.As  The changes of the number density of positive ions in the DBD plasma with pure oxygen at 1 kW RF power after the time of 1 μs and 1 ms are presented in Fig. 14.It clear that the positive ions have radially moved in the plasma from both dielectric and power electrode to the center of tube at the time of 1 ms although they were almost diffused in all the plasma region at 1 μs.
Moreover, the spatial variations of the number density of positive ions in the considered plasma discharge at different contents of argon in the gas mixture with the RF power of 1 kW are shown in Fig. 15.It can be observed that the number density of positive ions increases at the higher argon percentages in O 2 /Ar gas mixture.This  is owing to the higher cross section of momentum transfer for the atomic argon 56,57 , since collision energy of electrons decreases at the higher argon percentages.So, the recombination process between the electrons and positive ions will slightly occur in the plasma.It must be noted that, in addition to presence of O + and O + 2 , one of the most important positive ions is Ar + in this case of plasma (see Table S1) which its density will be increased by adding argon gas to the oxygen plasma and rising its contents.
In addition, the negative ions also have the same behavior as the positive ions (see Fig. 13) near the power electrode in spite of the fact that they show an opposite behavior in respect to the dielectric wall as seen in Fig. 16.As can be seen the negative ions are also tending toward the ground electrode similar to the electrons (see Fig. 10) while the positive ions are depleting around the ground electrode.All the charged species have depleted around the power electrode after 1 ms.In addition, the number density of negative ions is lower than the electron number density (see Fig. 10), both of them are lower than positive ions due to quasi-neutrality condition.
As presented in Fig. 16, the negative ions in the plasma are highly confined in the bulk region with a circular shape.The trapping effect of positive plasma potential is able to create this confinement of negative ions in the case of plasma.It should be mentioned that the electron attachment is mainly responsible for producing the negative species (see Table S1), and therefore the negative ions distribution has a similar behavior by the spatial distribution of electrons density in the considered plasma.The electronegativity of negative ions is equal to 0.14 by considering the maximum values of electrons (4.9 × 10 14 cm −3 ) and negative ions (0.7 × 10 14 cm −3 ) in pure oxygen plasma at 1 ms if the relation of α = n − /n e is used.
Meanwhile, the spatial changes of the number density of negative ions in the wire-to-multiwire DBD plasma at various contributions of argon in the O 2 /Ar gas mixture and 1 kW RF power is introduced in Fig. 17.It is obvious that the number density of negative ions decreases at the more argon contents in O 2 /Ar gas mixture.It must be noted that, due to the electronegativity of oxygen, the negative ions are created by the dissociative attachment and electron impact detachment reactions.However, this is an important loss mechanism for negative ions in all the oxygen plasma discharges.Although the negative ions density is affected by the plasma electrons density due to the dissociative attachment reactions, the amount of oxygen in the plasma discharge medium is reduced by increasing argon gas.Therefore, the number density of negative ions decreases at the more Ar contents in the O 2 /Ar gas mixture.
Despite of the fact that the negative ions and electrons move to the dielectric wall during the time, the positive ions are accumulated at the radial center of tube and go away from the power electrode and dielectric wall.Moreover, the spatial distributions of all species become slightly asymmetrical during the time as observed in Figs. 10, 14, and 16.This owing to the fact that the electron heating will be more local during the time which  results in the noticeable ionization near to the dielectric layer, and this causes the slightly asymmetrical distribution.Furthermore, the diffusion is reduced in this case of plasma because of working at atmospheric pressure.On the other hand, to reach a uniform and center-peaked distribution of charged species, the electron diffusion process and the loss of electrons in the attachment mechanism should be balanced in the oxygenate plasma that is not established in the wire-to-multiwire DBD plasma.
The profile of the electric field in the wire-to-multiwire dielectric barrier discharge at 1 kW RF power with pure oxygen and different argon added to the oxygen plasma is illustrated in Fig. 18.It can be seen that the electric field has a maximum value close to the power electrode and it gradually decreases with the distance from the power electrode.Additionally, the electric field is slightly increased at higher percentages of argon gas in the mixture due to the increase of charged particles in the plasma medium.

Conclusions
In this work, a new plasma jet was introduced by a wire-to-multiwire configuration and its operation in RF regime was numerically studied.Since it can assume that a wire consists of a large number of aligned points in a line, a wire is more singular than a surface.Therefore, the geometry used in designing this jet, i.e., wire-to-multiwire,

Figure 2 .
Figure 2. The simulation domain for wire-to-multiwire dielectric barrier discharge in the 2D model (the grid size is 2 nm).

Figure 4 .
Figure 4. Variation of the relative intensity of species in the discharge with the Ar concretion in the O 2 /Ar admixture.

Figure 5 .
Figure 5. Variations of (a) electron number density (with 4% error) and (b) rotation, vibration, and excitation temperatures as a function of Ar fraction in O 2 /Ar mixture at position 5 mm (with 4% error).

Figure 6 .
Figure 6.Temporal variations of electron (a) density and (b) temperature at different percentage of argon mixed by oxygen.

Figure 7 .Figure 8 .
Figure 7. Temporal variations of density profile for (a) O + and (b) O + 2 at different percentage of argon mixed by oxygen.

2 ( 48 )Figure 9 .
Figure 9. Temporal variations of density profile for (a) O and (b) O 3 at various percentages of oxygen mixed by argon.

Figure 11 .
Figure 11.Spatial variations of density profiles for electrons at different percentage of argon and 1 ms.

Figure 12 .
Figure 12.Spatial variations of temperature profiles for electrons at different percentage of argon and 1 ms.

Figure 13 .
Figure 13.Spatial variations of plasma (electron) frequency profiles for electrons at different percentage of argon at 1 ms.

Figure 14 .
Figure 14.Spatial variations of density profile for positive ions at 1 μs and 1 ms with pure oxygen.

Figure 15 .
Figure 15.Spatial variations of density profile for positive ions at different percentage of argon after 1 ms.

Figure 16 .
Figure 16.Spatial variations of density profile for negative ions at 1 μs and 1 ms with pure oxygen.

Figure 17 .
Figure 17.Spatial variations of density profile for negative ions at different percentage of argon after 1 ms.

Figure 18 .
Figure 18.Spatial variations of electric field profile at pure oxygen and different percentage of argon after 1 ms.

Table 1 .
The list of species considered in the model.

Table 2 .
Parameters used in the simulation.